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The effects of hydrogen impurities in the bulk and on the surface of aluminum are theoretically 
investigated. Within the framework of density functional theory, we have obtained the dependence 
on H concentration of the stacking fault energy, the cleavage energy, the Al/H surface energy and 
the Al/H/Al interface formation energy. The results indicate a strong dependence of the slip energy 
barrier in the [211] direction the cleavage energy in the [111] direction and the Al/H/Al interface 
formation energy, on H concentration and on tension. The dependence of the Al/H surface energy 
on H coverage is less pronounced, while the optimal H coverage is < 0.25 monolayer. The calcu- 
lated activation energy for diffusion between high symmetry sites in the bulk and on the surface is 
practically the same, 0.167 eV. From these results, we draw conclusions about the possible effect of 
H impurities on mechanical properties, and in particular on their role in embrittlement of Al. 



PACS numbers: 61.72.Qq, 61.72.Nn, 62.2.-X, 66.30.Fq, 68.43.Bc 
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I. INTRODUCTION 

In many technological applications of advanced mate- 
rials a crucial aspect of performance is the control of envi- 
ronmental effects, such as the presence of impurities. One 
such impurity is hydrogen, which pervades most metals 
and degrades their performance.!!! The interactions of H 
with lattice imperfections, such as dislocations, stacking 
faults, surfaces and microcracks, dominate its influence 
on the mechanical properties of a material. However, 
these interactions are far less well understood at a fun- 
damental level than the behavior of H in perfect crystals. 
Therefore, atomistic studies based on parameter-free, ab 
initio calculations are of great interest because they can 
provide accurate energetics for the various H-defect com- 
plexes and probe the microscopic physics responsible for 
the macroscopic behavior. The impurity-defect energet- 
ics are not only interesting by themselves, they can also 
be incorporated in more sophisticated models in order to 
make quantitative predictions for the macroscopic prop- 
erties of solids, in what has become known as multiscale 
simulations of materials. Bel 

The present study is motivated by the desire to shed 
light into H embrittlement of Al from an electronic struc- 
ture point of view. Experimentally, the presence of 
H in Al is associated with enhanced dislocation activ- 
ity which, ]j«ij]p|aps paradoxically, leads to a brittle rup- 
ture failure.tloO Theoretically, it has been shown recently 
within the framework of the Peierls-Nabarro model, that 
the presence of H in Al can dramatically enhaiice disloca- 
tion mobility and inhibit dislocation cross-slip.H However, 
the underlying atomic bonding features that give rise to 
such dislocation behavior have not been explored. In this 
paper, we show how H can change the nature of chemical 
bonding in Al, leading to. so-called hydrogen enhanced 
local plasticity (HELP).!i!'0 Moreover, we show that the 
cleavage energy, which represents the ultimate resistance 



to fracture, can be considerably reduced by H. 

Another important aspect of H behavior in Al is the 
thermodynamics of H in bulk Al and on its (111) surface, 
and the corresponding diffusion constants. The stability 
and mobility of H impurities in Al play an important role 
in HELP and H embrittlement of Al. For example, it is 
observed experimentally that HELP occurs only when 
the thermal diffusion of H in the lattice is fast enough 
to follow the motion of dislocations.!^! Accordingly, we 
calculate the diffusion energy barriers and the diffusion 
constants for H in the bulk and on the (111) surface of Al. 
We also examine the H diffusion process under uniform 
tensile strain to simulate the behavior at the dislocation 
core or near the crack tip region. 

The outline of this paper is as follows: we briefly 
describe our computational methodology in Sec. II. In 
Sec. HI we present the detailed results for the energetics 
of various defect structures involving H in the bulk and 
on the (111) surface of Al and attempt to understand 
some of the energetics from an electronic structure point 
of view. We discuss the physical consequences of our re- 
sults on the mechanical properties of Al and conclude in 
Sec. IV. 



II. METHODOLOGY 

The ah initio calculations we performed are 
based on denaty functional theory with the VASP 
implementatiorJS and ultra-soft pseudopotentials.l3 We 
have performed fc-point convergence studie.s in all cases 
using a uniform Monkhorst-Pack scheme.B From these 
studies we have determined that a grid consisting 
of 16x16x16 divisions in the Brillouin-Zone of the 
primitive unit cell of bulk fee Al, appropriately scaled 
for larger unit cells, is adequate for good convergence. 
The kinetic energy cutoff of 130 eV for pure Al yields 
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well converged results, whereas a higher cutoff of 350 
eV is needed in the presence of H atoms. We have 
also introduced a smearing of the Fermi surface by 
a temperature of 25 meV. With these computational 
parameters, the calculated lattice constant for bulk fee 
Al is a = 3.99 A and the bulk modulus is B = 83.2 GPa, 
determined by a Birjch-Murnaghan fit to the energy 
versus volume curve^ These values compare well with 
experimental valuesll^ of a = 4.05 A and B = 76.9 GPa, 
respectively. 

We next turn our attention to the atomic structures 
used to represent the physical systems of interest. In 
Fig. |l| we show the basic supercell which consists of 6 Al 
layers in the [111] direction. In all the calculations, we 
used the same number of layers in this direction to repre- 
sent the (111) surface or the interface between two semi- 
infinite slabs. We have used this supercell and multiples 
or distortions of it to calculate the generalized stacking 
fault energy curve and the surface and interface ener- 
gies. The unit cell with periodicity in the (111) plane 
equal to that of the bulk crystal will be referred to as 
the 1x1 cell. Multiples of the in-plane vectors, denoted 
as oi and 02 in Fig.l], were used to create larger super- 
cells for studying the effects of H concentration. We have 
used (qdi X 902) multiples of the basic cell with q = 1, 
and 2. In each supercell we included one H atom. 
These configurations correspond to H concentrations in 
the bulk of 14.3, 5.3 and 4.0 at.%, respectively. On the 
surface, they correspond to H monolayer (ML) coverages 
of e = 1.0, 0.333 and 0.25 ML, respectively. We also 
report a single surface calculation with a 4 x 4 unit cell, 
corresponding to H coverage 8 = 0.0625 ML, in order to 
establish the value of the H/Al surface energy in the low 
coverage limit. The H atom in each supercell was placed 
at the high symmetry interstitial sites, identified as the 
tetrahedral, 4-fold coordinated (T) or octahedral, 6-fold 
coordinated (O) position in the undistorted bulk config- 
uration, or as the fee (F), hep {H) and atop {A) sites 
on the (111) surface, all shown in Fig. |^. Calculations 
with the H atom in positions between the high symmetry 
sites in the bulk and on the surface were used to obtain 
diffusion energy barriers. 

Distorting the q x g supercell by increasing the hor- 
izontal or vertical components of 03, the lattice vector 
which in the undistorted case lies along the [111] direc- 
tion, produces configurations that generate the general- 
ized stacking fault energy or the cleavage energy. These 
two distortions are referred to as d and h respectively, and 
are given in their natural units of a/\/6 and a/\/3. In 
these units, d = 1 corresponds to the intrinsic stacking 
fault configuration, d = 2 corresponds to the so-called 
run-on configuration (in which two Al atoms on either 
side of the slip plane are exactly above and below each 
other) , and c? = 3 corresponds to another ideal configura- 
tion identical to d = 0. Similarly, h = 1 corresponds to a 
separation between the two slabs equivalent to a missing 
(111) layer. 

For each of these configurations, all atoms except the 
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FIG. 1: Schematic representation of the geometry used in the 
total-energy calculations. The basic supercell is shown in side 
and top views. The dashed lines outline the undistorted su- 
percell. It consists of 6 layers in the [111] direction and has 
repeat vectors in the (111) plane equal to the ideal crystal 
primitive lattice vectors, denoted by a\,a2. The third vector, 
as, is along the [111] direction in the undistorted cell. Dis- 
tortions of the cell by adding to as components in the [211] 
direction (denoted by d, given in units of a/\/6), or in the 
[111] direction (denoted by h, given in units of a/-\/3), lead to 
configurations relevant to the generalized stacking fault en- 
ergy surface or to the cleavage energy; such a distortion for 
d = 1,/i = 1 is illustrated. The large white, grey and black 
circles indicate the positions of the Al atoms; all Al atoms 
contained in a unit cell are shown in the side view, but only 
selected planes of atoms are shown in the top view. The 
smaller black circle indicates the high symmetry positions of 
the H atom in the bulk (tetrahedral and octahedral), and on 
the (111) surface (F for the fee site, H for the hep site and A 
for the atop site). 



innermost two layers of the Al slab were fully relaxed 
via the conjugate-gradient method, so that the magni- 
tude of the calculated forces on the atoms was less than 
0.03 eV/A. For the calculations of the energy barriers for 
diffusion, the coordinates of the H atom are held fixed, 
either in all three directions for bulk diffusion, or in the 
lateral surface directions for surface diffusion. For the 
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FIG. 2: The generalized stacking fault energy 7gsf for Al 
with H impurities at zero tension opening, h — 0, as a function 
of the slip d in the [211] direction for the 1x1 supercell, cor- 
responding to H concentration of 14.3 at.%: circles represent 
the energy for the tetrahedral site, squares for the octahedral 
site (lines are fits intended as guide to the eye) . 

calculation of the generalized stacking fault energies us- 
ing the distorted bulk supercells (see below), the H atom 
was placed initially close to the interpolated tetrahedral 
or octahedral positions and allowed to relax to the near- 
est local energy minimum. We report energy differences 
between various configurations in eV and surface energies 
in Jm~^, in order to comply with conventions in the lit- 
erature and make our results easily comparable to other 
published work. 



III. ENERGETICS OF H IMPURITIES IN AL 

A. Generalized stacking fault energies 

The generalized stacking fault (GSF) energy, denoted 
by JGSF, is defined as the energy cost per unit area for 
sliding two semi-infinite slabs relative to each other along 
a particular plane by a certain vector d. The energy sur- 
face generated by spanning the allowed values of d con- 
tains several important features relevant to the mechan- 
ical properties of solids and in particular to their brittle 
versus ductile behavior. For an fee metal like Al, the 
most interesting portion of the GSF energy surface is the 
path along the [211] direction on the (111) plane. This 
path includes both the intrinsic stacking fault as well as 
the unstable stacking fault, corresponding to d = 1 and 
d = 0.6, respectively. The intrinsic stacking fault en- 
ergy, denoted as 7^/, along with the elastic properties 
of the material, determiniL the separation distance be- 
tween partial dislopationsllJ which controls the mobility 
of the dislocations. EJ The unstable stacking fault energy, 
denoted as "/us , represents the energy barrier for disloca- 
tion nucleation from a crack tip, which is related to the 



tendency for brittle or ductile behavior of the material.113 
The val|Ues-pf jgsf for pure Al have been published 
elsewhere. E3lJ In the present work we have repeated 
these calculations to obtain a consistent set of numbers 
with the computational parameters and methodology 
adopted here. The values of the important configurations 
obtained by the present calculations are: 7„s = 0.182 
Jm~^, and 7s/ — 0.134 Jm~^, both for h — 0. Including 
tension opening {h 7^ 0), reduces these values dramati- 
cally. For example, for /i = 0.1, we obtained 7„s = 0.094 
Jm-^ a 50% reduction, and 7^/ = 0.092 Jm'^, a 30% 
reduction. 

The values of "/gs f in the presence of H depend on the 
position of the H atom in the lattice and the H con- 
centration. We next examine these two contributions 
separately. We consider first our findings for the high- 
est H concentration, 14.3 at.%, which corresponds to 
one H atom in a 1x1 supercell. The results, shown in 
Fig. 1^, indicate that there is a cross-over at approximately 
d = 0.5 in site preference for the H atom. This cross-over 
in site preference significantly reduces the unstable stack- 
ing energies to jus = 0.097 Jm~^, a ^ 50% reduction, and 
the intrinsic stacking fault energy to 73/ = 0.073 Jm~^, 
also a ~ 50% reduction. The reason for such energy re- 
duction is that the volume available for the interstitial 
H atom situated at the original octahedral site decreases 
during the slip, while it increases at the tetrahedral site. 
The effects of tension on 7gs_f in the presence of H were 
also calculated and found to be very similar to those for 
pure Al, as far as the relative energy decrease is con- 
cerned. 

To study the dependence of 7^^ and jsf on H concen- 
tration, we have computed these energies for several H 
concentrations. The results are given in Table I. The gen- 
eral trend for both energies is to increase with decreasing 
H concentration. An exception to the trend is the highest 
H concentration at 14.3 at.% (from the 1x1 supercell), 
which has higher values for 7„s and 7^/ than the next 
lower H concentration of 5.3 at.% (from the x \^ su- 
percell) . We believe that this has to do with the fact that 
in the 1x1 supercell at fixed volume, the optimal ionic 
bonding distances between H and Al atoms cannot be 
satisfied, and therefore the system cannot attain a struc- 
ture with a reasonably low energy. To investigate this 



TABLE I: The unstable stacking energy, 7„s, stacking fault 
energy, 7^/, and cleavage energy, 'yd, for the H/Al system as a 
function of H concentration without volume relaxation. The 
ratio 7c;/7us is also included. 
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1x1 


14.3 


0.097 


0.073 


0.930 


9.6 




5.3 


0.089 


0.071 


1.611 


18.1 


2x2 


4.0 


0.136 


0.074 


1.680 


15.6 


1x1 


0.0 


0.182 


0.134 


1.934 


10.6 
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FIG. 3: The bonding valence charge density on the(l21) plane 
for pure Al (left) and Al+H (right) systems at d = 0. The 
horizontal direction is [101] and the vertical direction is [111]. 
The fictitious slip plane is shown as a dashed line and the 
positions of atoms are indicated by the corresponding labels. 
The arrow indicates the buildup of valence charge density 
corresponding to covalent bonding across the slip plane. 



conjecture, we have also examined thoroughly the effect 
of volume relaxation in the 1x1 supercell. We find that in 
this supercell there is actually a low energy configuration 
at the run-on position, d = 2, which is lower in energy by 
0.615 Jm^^ than the undistorted configuration of = 0, 
and involves an increase in the volume by 18%. In this 
configuration, the H atom lies exactly between the two 
Al atoms on either side of the slip plane, forming strong 
ionic bonds across the interface, which compensates the 
energy loss due to the distortion of the Al lattice. 



B. Cleavage energy 

Another important energy for the mechanical behav- 
ior of a solid is the cleavage energy, 7ci, defined as the 
energy cost per unit area to separate the solid into two 
semi-infinite halves by creating two surfaces. For the 
pure Al case, 7c/ = 27^, with 7^ the energy of the newly 
created surface. We have calculated the cleavage energy 
with the various supercells as the vertical component of 
03 is increased up to ft, = 4, corresponding to a separation 
of about 9 A. The H atom was placed at the energetically 
preferred tetrahedral site. With the introduction of H to 
the system, the cleavage energy dramatically decreases 
by as much as 50% at H concentration of 14.3 at.% (see 
Table I) . The decrease of the cleavage energy is approxi- 



mately proportional to the H concentration. 

In order to elucidate the origin of the reduction in the 
GSF energy and the cleavage energy in the presence of 
H, we examine the bonding charge density on the (121) 
plane for pure Al and for Al-t-H (14.3 at.%) at d = 0, 
which is shown in Fig. ^. The bonding charge density is 
defined as the difference between the valence charge den- 
sity in the solid and the superposition of neutral atomic 
valence charge densities placed at the lattice sites. The 
positive (negative) bonding charge density represents the 
net gain (loss) of charge as the atoms are brought to- 
gether to form the solid. The contour graph is shaded in 
such a way that regions with higher value of charge den- 
sity are lighter. By examining the bonding charge density 
of Al and Al-I-H, we find extended covalent bonding (in- 
dicated by the arrow) in Al across the slip plane, which 
is dramatically weakened in the presence of H. In fact, 
the H atom depletes the Al bonding charge from the in- 
terstitial region and the regions across the slip plane to 
form ionic bonding between the H sites and the nearest 
Al sites above it. As a consequence, the cohesive strength 
across the slip plane is reduced by the presence of H, giv- 
ing rise to the lower cleavage energies. More importantly, 
since the strength of the ionic bonding between the pos- 
itively charged Al plane and negatively charged H plane 
is not sensitive to the relative sliding between the two 
planes, the sliding energy barrier is greatly reduced, and 
the GSF energy surface becomes much smoothcifl in the 
presence of H. This is contrasted to the pure Al case, 
where the covalent bonding among Al atoms across the 
slip plane is very sensitive to the local bonding distor- 
tions and consequently thejJjSF energy is higher and has 
more pronounced features. lj Although these calculations 
concern Al, we believe that the results are also applicable 
to other metals whose electronegativity is lower than H. 

In Table I we also give the ratio of 7ci to 7,1^ for the 
various H concentrations with pure Al as the reference 
point. The value of this ratio is indicative of the tendency, 
of the material to exhibit brittle or ductile behavior.EBl 
From this comparison we infer that at modest H con- 
centrations the system has increased ductility, which is 
consistjfint with the experimental observations concerning 
HELPjd but at the highest H concentration considered 
the system may become less ductile. The anomaly of the 
7c/ llus ratio at the highest H concentration is related to 
the anomalous behavior of at this concentration, as 
noted earlier. 



C. Absorption and diffusion of H in bulk Al 

In view of the importance of the thermodynamics of H 
in bulk Al and the kinetics of H transport in the presence 
of defects, we have performed additional calculations for 
the energetics of H absorption and the H diffusion energy 
barrier in bulk Al. For these calculations we employed 
a 32-atom supercell of the bulk crystal, which is a mul- 
tiple of the conventional simple cubic cell by a factor of 
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Tetrahedral Octahedral Tetrahedral 

H position 

FIG. 4: The relative energy for motion of a H atom from the 
tetrahedral position (0) to the octahedral position (1) in the 
32-atom bulk supercell. Circles represent the energy without 
tensile strain, and squares correspond the energy with 5% 
tensile strain. 

2 in each direction (hence, we refer to it as the 2x2x2 
supercell). We have investigated the two high symmetry 
interstitial sites of a single H impurity, which corresponds 
to a H concentration of 3.03 at.%. The absorption en- 
ergy, Eab, was obtained with reference to the cohesive 
energy of crystalline fee Al, Ec{A\), and using a gas of 
H2 molecules, whose binding energy is £^b(H2), as a reser- 
voir for the H atoms: 

Eab = E.in/Al) - N,,E,iAl) - ^Eb{n2), (1) 

where £'c(H/Al) is the calculated cohesive energy of the 
Al supercell configuration with one H impurity and N^c 
is the number of Al atoms in the supercell {Nsc = 32 in 
the present case). The binding energy of the H2 molecule 
was calculated to be —6.697 eV in vacuum, using a cubic 
cell with side equal to 24 A, and the same computational 
parameters as for the H/Al system. We find -E^f, — 

-0.222 eV for the tetrahedral site and = -0.152 

eV for the octahedral site. These results show that the 
incorporation of H in bulk Al, starting with an Al crystal 
and H2 gas, is a thermodynamically exothermic process. 
Therefore the H impurity is thermodynamically stable in 
bulk Al. 

Having established the stability of H in bulk Al, we in- 
vestigated the mobility of H in the Al lattice. As alluded 
to earlier, the diffusion rate of H in the lattice determines 
whether H-enhanced local plasticity can occur or not. 
HELP takes place only when the diffusion of H atoms is 
fast enough to allow them to redistribute around the core 
of a moving dislocation (dynamic trapping! and thereby 
continuously minimize the system energyEI H diffusion 
can also affect the kinetics of crack propagation, the 
strain rate dependenoe-pf H embrittlement, and the rate 
of hydride formationJi3 Since H jumps between nearby 



interstitial sites in bulk Al, we calculated the activation 
energy of H diffusion between the closest tetrahedral and 
octahedral sites. In Fig. ^ we show the energy as a func- 
tion of relative position of the H atom moving between 
the tetrahedral (0 and 2) and the octahedral site (1). 
Another special site is the mid-point between the tetra- 
hedral and the octahedral sites. We define the energy 
of the mid-point and of the octahedral site relative to 
the tetrahedral site as ei and £2, respectively. The diffu- 
sion energy barrier for unstrained bulk Al is Sb — 0.162 
eV, which compares well with the experimental valuec£l 
of 0.168 eV. The energy barrier configuration happens 
to coincide with the mid-point between the octahedral 
and tetrahedral sites in the unstrained crystal, that is 
Eb = Si. Since we are also interested in the interaction of 
H with lattice imperfections, such as dislocations, micro- 
cracks, etc., we considered how H diffusion is affected by 
the presence of such defects. A simple way to simulate 
this effect is to apply strain to the system. In this study 
we concentrate on how tensile strain affects H diffusion 
in Al, since this is the type of strain field usually found 
around an edge dislocation or a crack tip under model I 
loading, which are both relevant to H embrittlement of 
Al. For the range of tensile strains studied, the tetra- 
hedral site has the lowest energy. When the hydrostatic 
tensile strain is small (< 3%), the energy barrier is lo- 
cated at the mid-point between the tetrahedral and the 
octahedral sites, i.e., Sb = £i- Moreover, we find that 
for small strain the diffusion energy barrier drops mono- 
tonically as the strain increases, and is reduced to 0.150 
eV for 3% expansion. This result is important because 
it shows that not only H prefers to stay in slightly en- 
larged interstitial regions such as dislocation cores and 
crack tips, but that it can also move more easily within 
such regions. On the other hand, if the tensile strain is 
large (> 3%), the octahedral site becomes energetically 
unstable and represents the energy maximum where the 
energy barrier for diffusion is located, that is = 62- 
One example for such diffusion energy profile is shown 
in Fig. m corresponding to 5% strain. The instability 
of the octahedral site arises from the unfavorable bond 
length between H and Al atoms, 2.1 A, which is much 
larger than the preferred ionic bond length of about 1.8 
A. The values of ei and £2 as a function of tensile strain 
are summarized in Table II. Noticing that the tetrahe- 
dral site is always energetically favorable regardless of the 
strain, we also calculated the relative energy of H at the 
mid-point of the direct line between two adjacent tetra- 
hedral sites defined as £3 as a function of tensile strain. 
We find this energy difference to be always higher than 
the corresponding values of £1 and £2. Therefore we have 
confirmed that the tetrahedral-octahedral-tetrahedral se- 
quence is the preferred diffusion path for H in bulk Al. 
The various values of £3 as a function of tensile strain are 
also listed in Table II. 

The calculated energy profiles allow us to estimate the 
bulk diffusion constant, 

Db = i^bll exp[-eb/kBT] (2) 
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where Vb and It are the attempt frequency and hopping 
length for bulk diffusion. Approximating the energy dif- 
ferences near the equilibrium tetrahedral configuration 
by a second order polynomial in the distance, we find an 
attempt frequency i^b — 0.8 x 10^^ sec~^, while the hop- 
ping length between equivalent sites is lb = 0.948 a. The 
diffusion constant at room temperature (300 K) is esti- 
mated to be 1.78 X 10^^^ m^s^^ in the unstrained crystal. 
Assuming that the values of the attempt frequency and 
hopping length are not significantly affected by strain, 
the value of the bulk diffusion constant at room temper- 
ature and for 3% and 5% tensile strain is 3.46 x 10^^^ 
m^s~^ and 1.31 x 10"-'^^ m^s~^, respectively. 



D. Adsorption and diffusion of H on AI(lll) 

For H embrittlement of metals, it has been experimen- 
tally observed that the fracture surface is-along the slip 
plane, where shear localization occurs.ElEZi Apparently, 
adsorption and diffusion of H on the fresh fracture surface 
play an important role on the kinetics of crack propaga- 
tion and the embrittling effect of H. The critical ener- 
getics that are relevant to the adsorption and diffusion 
process not only can provide insight into the problem, 
they can also be used in an empirical analysis which cau 
deal with the macroscopic aspects of this phenomenon.^ 

For this reason, we have considered the adsorption of a 
H atom on the (111) surface of Al, in configurations cor- 
responding to 1x1, -s/S X ^/S, and 2x2 surface unit cells, 
or H coverages in the range 9 e [0.25, 1.0] ML. The H 
atom was placed at the high-symmetry F, H and A points 
(see Fig. |^), as well as points at regular intervals between 
them, to determine the lowest energy configuration and 
the energy barrier for surface diffusion. By analogy to 
the definition of Eq. (0) , the adsorption energy Bad is de- 
fined as the energy of the configuration with a H atom on 
the (111) Al surface, relative to the same surface with- 
out H and using a gas of H2 molecules as a reservoir 
for the H atoms. For a full monolayer of H (O = 1.0) 



TABLE II: The energy of special points for H diffusion in bulk 
Al, relative to the tetrahedral position: ei is the energy of the 
mid-point between the tetrahedral and octahedral positions, 
£2 is the energy of the octahedral position and £3 is the en- 
ergy of the mid-point between adjacent tetrahedral sites. The 
asterisks denote the energy barrier for diffusion. 



Strain 


£1 


£2 


£3 


(%) 


(cV) 


(eV) 


(eV) 





0.167 * 


0.071 


0.359 


1 


0.161 * 


0.094 


0.332 


2 


0.156 * 


0.116 


0.306 


3 


0.150 * 


0.137 


0.280 


4 


0.143 


0.157 * 


0.255 


5 


0.134 


0.175 * 


0.229 



on the (111) Al surface, corresponding to the 1x1 cell, 
the adsorption energy for the various positions along the 
[211] direction is given in Fig. |. The fee site (F) is the 
energetically preferred position with Bad — —0.085 eV, 
whereas the hep site (H) has nearly zero adsorption en- 
ergy, and the atop site (A) is energetically unfavorable. 
We find that when displacing the H atom from the F 
toward the H and A positions, the Al atoms near the 
surface are also displaced in the [211] direction, so as to 
maintain the high coordination of the H atom to the ex- 
tent possible. This is especially pronounced at the bridge 
position, which is half way between the F and the H sites 
(see Fig. |). 

From the results of these calculations we conclude that 
the diffusion of H on the (111) Al surface will follow a 
zig-zag path between successive F and H sites. We find 
that the energy difference between the H and F sites 
is a reasonable approximation for the diffusion energy 
barrier, e^, within the numerical uncertainty inherent in 
the calculations. With this in mind, we have calculated 
the diffusion energy barrier as a function of H coverage, 
using multiples of the 1x1 surface unit cell which is 
given in Table III. It is clear from this Table that the 
energy barrier for surface diffusion of an isolated H atom 
is approximately 0.163 eV, practically the same as that 
for bulk diffusion. We can also obtain an estimate of the 
surface diffusion constant using the energy as a function 
of H position from the calculation of the 1x1 unit cell: 

Ds = iyslUM~£s/kBT] (3) 

where i>s and Is are the attempt frequency and hopping 
length for surface diffusion. Using the same procedure 
as for bulk diffusion, we find an attempt frequency Vg = 
0.7 X 10^'^ sec~^, the hopping length between equivalent 
sites is Is = 0.707 a, and the diffusion constant at room 
temperature (300 K) is 1.02 x 10"^^ m^s"^ 

Finally, we note that the presence of H on the (111) 
Al surface reduces the surface energy considerably. In 
order to quantify this observation, we report in Table 
III the calculated adsorption energy. Bad, for H atoms 
at various coverages. The reduction in surface energy in 
the presence of H, A'js, is obtained by converting the 
adsorption energy to a surface energy and subtracting 
from it the corresponding surface energy of pure Al. We 
find that this reduction in surface energy is a function 
of coverage, and it increases as the coverage decreases 
(Table III). Finally we have performed one additional 
calculation in a larger 4x4 surface unit cell in order to 
determine to what extent this trend continues for lower 
coverages; the result is included in Table III. We conclude 
that the H coverage which gives the largest reduction in 
surface energy is in the range 0.0625 < & < 0.25 ML. 

It should be pointed out that A'js is a very impor- 
tant material parameter in determining the tendency of 
impurity-induced intergranular fracture. More specifi- 
cally, according to the thermodynamic theory developed 
by Rice and Wangjla the potency of a segregating im- 
purity in reducing the Griffith work of a brittle grain 
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0.15 




H position 

FIG. 5: The adsorption energy for H on the (111) Al surface 
within the 1x1 surface unit cell corresponding to H coverage 
O — 1 ML, for different positions in the [211] direction: F — 
fee site, H — hep site, A = atop site. 



boundary separation is a linear function of the difference 
Ajgb — Ajs, that is, the difference between the segrega- 
tion energy of that impurity at a grain boundary and at 
a free surface. A smaller reduction in surface energy (a 
less negative Ajs) indicates a weaker tendency for brittle 
intergranular fracture. Based on the fact that our calcu- 
lated Ajs is more than an order magnitude-aaaller than 
typical values for intergranular fracture we infer 
H-induced fracture in Al to be of a transgranular nature. 
Of course, the definite determination of the tendency will 
also depend on Ajgb, which is not available to us. 

It is also instructive to compare this reduction in sur- 
face energy to the energy required to form Al/H/Al inter- 
face structures corresponding to the bulk unit cells with 
one H impurity discussed earlier. In Table III we give the 
absorption energy Eab obtained from the corresponding 
unit cells employed in the bulk calculations. As before, 
we use the pure Al bulk as a reference system and a gas 
of H2 molecules as the reservoir for H atoms. The con- 
figurations used for these calculations correspond to the 
formation of a planar interface between two (111) planes 
of bulk Al separated by an ordered layer of H atoms. We 
have calculated the interface formation energy, Aji^t , as 
a function of H content at the interface, expressed in ML 
of H; the results are given in Table III. Consistent with 
the calculations we presented earlier, this formation en- 
ergy is positive for large H concentrations (recall the large 
positive absorption energy for a H atom in the 1x1 bulk 
supercell). However, with decreasing H concentration at 
the interface we expect that this formation energy will be 
reduced, and for small enough concentrations it should 
be negative, corresponding to the absorption energy for 
an isolated H atom in bulk Al, which we found to be 
—0.222 eV. Indeed, Ajint becomes negative for the 2x2 
supercell. Note that in this supercell the shortest dis- 
tance between H atoms on the (111) plane is aV^, which 



is shorter than the distance between the H impurities in 
the 32-atom bulk supercell, equal to 2a. 



IV. DISCUSSION AND CONCLUSIONS 

The following picture emerges from the calculations 
reported above. In a system consisting of a crystalline 
Al phase in equilibrium with a II2 gas, H atoms will be 
adsorbed on the (111) Al surfaces, the natural cleavage 
planes, in order to lower the surface energy. The equi- 
librium coverage will be in the range 0.0625 < O < 0.25 
ML. Diffusion of H on the surface and in the bulk is rel- 
atively fast. Using the values for the diffusion constants 
we determined earlier, we find that at room temperature 
the length scale for surface diffusion in a time interval of 
1 sec is ~ 0.3 /im, while for bulk diffusion it is ~ 7 /im. 
Over such distances the equilibrium surface H coverage 
and bulk H content will be dictated by thermodynamic 
considerations. The calculated bulk absorption energy 
for the tetrahedral site, -E^P = —0.222 eV, indicates 
that energetically it is possible for H to end up in the 
bulk. However, this is not the preferred configuration. 
In fact, the calculated surface adsorption energies per H 
atom are lower than the bulk absorption energies. For 
example, from the calculation of the 4x4 surface unit 
cell with one H atom at the F position, we find a surface 
adsorption energy E^^^^ = —1.907 eV, significantly lower 
than the bulk absorption energy. In this configuration 
the H atom can be considered as an isolated atom on 
the Al surface. It appears from these calculations that 
if there is any H in bulk Al and if the system is allowed 
to equilibrate with the surface, H tends to diffuse out 
and remains on the surface. This conclusion rests on the 
assumption that a clean and atomically flat (111) Al sur- 



TABLE III: The energy difference between the two high sym- 
metry positions F and of a H atom on the Al(lll) surface, 
identified as the surface diffusion activation energy Es ; the ad- 
sorption energy E^^ of H atoms at the energetically preferred 
fee (F) site on the (111) surface and the corresponding reduc- 
tion in surface energy A'js; the absorption energy -E^P of H 
atoms at the energetically preferred tetrahedral (T) site in 
the interface between two (111) planes and the corresponding 
interface formation energy A^int- All quantities are given as 
functions of H coverage O in monolayers (ML). The last line 
gives the corresponding results for the 32-atom bulk supercell 
of the conventional cubic cell. 



Supercell 


e 


es 




A7, 




Ajint 




(ML) 


(eV) 


(eV) 


(Jm-2) 


(eV) 


(Jm-2) 


1x1 


1.0 


0.092 


-0.085 


-0.138 


+0.383 


+0.586 




0.333 


0.157 


-0.365 


-0.187 


+0.067 


+0.034 


2x2 


0.25 


0.163 


-0.489 


-0.188 


-0.039 


-0.015 


4x4 


0.0625 




-1.907 


-0.183 






2x2x2 




0.167 






-0.222 
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face is available, which is usually not the case in reality 
due to the tendency of Al surfaces to oxidize. The pres- 
ence of an oxide on the surface will completely alter the 
thermodynamic balance, making it possible for H atoms 
to remain in the bulk. Moreover, we have not investi- 
gated the possibility of surface vacancies or other defects 
on the surface, which can also change the thermodynamic 
picture. The presence of such defects can only lower the 
Al surface energy, while the binding of H atoms at such 
defects may not be preferred over binding on the flat sur- 
face, as the values of and E^^'' reported in Table III 
suggest. Thus, the presence of defects on the Al surface 
may also suffice to make the incorporation of H atoms in 
the bulk thermodynamically stable. 

On the other hand, our calculations for H in bulk Al 
show that H atoms feel an effective repulsive interaction 
for distances shorter than twice the primitive lattice con- 
stant, a/\/2. This is evident from the Al/H/Al interface 
energies Aji„t, reported in Table III, which are all posi- 
tive except for the 2x2 supercell. Therefore, H atoms in 
bulk Al cannot form dense clusters but have to be apart 
from each other by a distance at least a/\/2. This result 
is significant because it spells out the importance of H-H 
interactions and casts doubt on studies that ignore such 
interactions. The energetically favorable configuration 
which we found for high H concentrations in the bulk, 
with a H atom between two Al atoms directly above and 
below it in the [111] direction, also has interesting im- 
plications. Since this structure has a lower energy than 
the undistorted crystal configuration with equal H con- 
centration but also involves a large volume relaxation 
(18%), we conclude that if there are voids or other de- 
fects which give rise to tensile strain, H atoms will be 
preferentially bound to those sites, such as the cores of 
edge dislocations. In fact, a recent study showed that the 
binding energy of H to the core of an edge dislocation is 
much larger than that of a screw dislocation,0 which is 
in line with our observation here. These results have an 
important consequence for dislocation motion: the edge 
dislocation needs to turn into a screw dislocation in or- 
der to cross-slip, a process that will be hindered by the 
binding of H atoms to the edse. dislocation, consistent 
with experimental observations.Eil More importantly, this 
H-inhibited cross-slip will give rise to slip planarity and 
possibly, shear localization, the two most important el- 
ements to understand H embrittlement in terms of the 
HELP mechanisniE 

Finally, we turn our attention to the effects of H on 
the intrinsic stacking fault energy 7^/, unstable stacking 
fault energy jus, and cleavage energy, jci- As already 
noted, the presence of H reduces all these quantities rel- 
ative to their values in pure Al. The ratio jci/"fus as a 
function of impurity content has baen employed to dis- 
cern brittle versus ductile responsccS The simple physical 
picture behind this argument is that a low value of this 
ratio indicates a preference for cleavage rather than dis- 
location generation at a crack tip which is controlled by 
the value of 7us;li3 this behavior is associated with brit- 



tle failure. Conversely, a high value of this ratio indicates 
the preference for dislocation generation at a crack tip, 
a behavior associated with ductile response. While this 
picture may be overly simplistic for quantitative anal- 
ysis, it does give insight about general trends and can 
even lead to useful predictions that have been verified 
experimentally.c3 For H in bulk Al, this ratio is higher 
than the value in pure Al except at the highest H concen- 
tration considered here (14.3 at %, corresponding to the 
1x1 supercell) , clearly suggesting that H embrittlement 
in Al takesriplace as plastic rupture rather than as brittle 
separation.EJ In fact, recent worka indicates that H in Al 
can indeed lead to etihanced local plasticity, a precursor 
to H embrittlement .Id 

The lower value of the intrinsic stacking fault energy 
in the presence of H suggests a larger separation of the 
partial dislocations in Al, which could hinder the dislo- 
cation cross-slip since the partial dislocations must be 
constricted before cross-slip can take place. But a mora 
detailed analysis based on the Peierls-Nabarro modelEI 
shows that, even though the stacking fault energy is low- 
ered by the presence of H, the partial dislocations are not 
split any further than in pure Al, while the core width 
of the dislocations is increased significantly giving rise to 
enhanced dislocation mobility. 

To summarize, we have performed density functional 
theory calculations to study the energetics of H impu- 
rities in bulk and on the (111) surface of Al. We have 
obtained the dependence of the stacking fault energy and 
the cleavage energy, as well as the Al/H surface energy 
and the Al/H/Al interface formation energy, on H con- 
centration. The results indicate that there is a strong 
dependence of the GSF energy in the [211] direction, the 
cleavage energy in the [111] direction and the Al/H/Al 
interface formation energy, on H concentration and on 
tensile strain. We are able to explain the H-induced 
reduction of the stacking fault energy and cleavage en- 
ergy in Al from an electronic structure point of view, 
and conjecture that such reduction can also take place 
in other H-metal systems. It is found that the depen- 
dence of the Al/H surface energy on the H coverage is 
less pronounced, with the optimal coverage being < 0.25 
monolayer. The calculated activation energy for diffu- 
sion between high symmetry sites in the bulk and on the 
surface is practically the same, 0.167 eV, in good agree- 
ment with experimental measurements. Although our 
calculations reported here provide strong theoretical ev- 
idence for the HELP mechanism, they are not able to 
answer how HELP eventually leads to H embrittlement. 
Nevertheless, we believe that our work sets the stage for 
developing a comprehensive theory of H embrittlement, 
which most likely will necessitate a multiscale framework. 
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